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Abstract 

Background: Bitter acids {e.g. humulone) are prenylated polyketides synthesized in lupulin glands of the hop plant 
{Humulus lupulus) which are important contributors to the bitter flavour and stability of beer. Bitter acids are fornned 
from acyl-CoA precursors derived from branched-chain amino acid (BCAA) degradation and C5 prenyl diphosphates 
from the methyl-D-erythritol 4-phosphate (MEP) pathway. We used RNA sequencing (RNA-seq) to obtain the 
transcriptomes of isolated lupulin glands, cones with glands removed and leaves from high a-acid hop cultivars, 
and analyzed these datasets for genes involved in bitter acid biosynthesis including the supply of major precursors. 
We also measured the levels of BCAAs, acyl-CoA intermediates, and bitter acids in glands, cones and leaves. 

Results: Transcripts encoding all the enzymes of BCAA metabolism were significantly more abundant in lupulin 
glands, indicating that BCAA biosynthesis and subsequent degradation occurs in these specialized cells. 
Branched-chain acyl-CoAs and bitter acids were present at higher levels in glands compared with leaves and cones. 
RNA-seq analysis showed the gland-specific expression of the MEP pathway, enzymes of sucrose degradation and 
several transcription factors that may regulate bitter acid biosynthesis in glands. Two branched-chain 
aminotransferase (BCAT) enzymes, HIBCATl and HIBCAT2, were abundant, with gene expression quantification by 
RNA-seq and qRT-PCR indicating that HIBCATl was specific to glands while HIBCAT2 was present in glands, cones 
and leaves. Recombinant HIBCATl and HIBCAT2 catalyzed forward (biosynthetic) and reverse (catabolic) reactions 
with similar kinetic parameters. HIBCATl is targeted to mitochondria where it likely plays a role in BCAA catabolism. 
HIBCAT2 is a plastidial enzyme likely involved in BCAA biosynthesis. Phylogenetic analysis of the hop BCATs and 
those from other plants showed that they group into distinct biosynthetic (plastidial) and catabolic (mitochondrial) 
clades. 

Conclusions: Our analysis of the hop transcriptome significantly expands the genomic resources available for this 
agriculturally-important crop. This study provides evidence for the lupulin gland-specific biosynthesis of BCAAs and 
prenyl diphosphates to provide precursors for the production of bitter acids. The biosynthetic pathway leading to 
BCAAs in lupulin glands involves the plastidial enzyme, HIBCAT2. The mitochondrial enzyme HIBCATl degrades 
BCAAs as the first step in the catabolic pathway leading to branched chain-acyl-CoAs. 



Background 

The female inflorescences (cones') of Humulus lupulus L. 
(hop, Cannabaceae) contain prenylated acylphlorogluci- 
nols and prenylchalcones (e.g. xanthohumol) that are 
important for the brewing industry and have potential 
medicinal uses [1-3]. The acylphloroglucinols, referred to 
as bitter acids because they contribute the characteristic 
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bitter flavour of beer, include the a-acid humulone and its 
acyl-side chain variants cohumulone and adhumulone 
(Figure la). Antimicrobial p-acids (e.g. lupulone, colupu- 
lone and adlupulone), which differ from a-acids by the 
presence of an additional dimethylallyl prenyl group, are 
also abundant. Hop cultivars vary widely in their content 
and composition of bitter acids with some "super- alpha" 
cultivars containing greater than 20% humulone by dry 
weight in cones [4]. 

The bitter acids and other prenylated polyketides 
occur mainly in lupulin glands, which are glandular 
trichomes present on hop cones and leaves [5]. 
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Figure 1 Biosynthesis of iiop bitter acids from branched-chain amino acid (BCAA) precursors, a) The acyl side-chains of humulone/ 
lupulone are derived from leucine, cohumulone/colupulone from valine and adhumulone/adlupulone from isoleucine via the corresponding 
acyl-CoA thioesters. The carbon chain of the BCAAs is shown in red to indicate their incorporation into the a- and |3-acids. b) Schematic 
depiction of the major pathways contributing to the biosynthesis of the major a-acid product, humulone, in lupulin glands. Pathways are 
coloured to match Figure 2b: BCAA biosynthesis, purple; BCKDH complex, brown; MEP pathway, orange; malonyl-CoA biosynthesis, blue; bitter 
acid biosynthesis, green. The movement of metabolic intermediates between cellular compartments is indicated by dashed line. Note: the 
subcellular localization of the putative humulone synthase enzyme is not known. 
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Glandular trichomes typically have little or no photosyn- 
thetic activity and often contain only a few highly active 
biosynthetic pathways, which has made transcriptome 
analysis of isolated trichome secretory cells useful for 
the identification of the genes involved in the produc- 
tion of specialized metabolites [6,7]. Enzymes involved 
in the biosynthesis of xanthohumol and terpenoids in 
hop lupulin glands [5,8] and cannabinoids in glandular 
trichomes of the closely related species Cannabis sativa 
(Cannabaceae) [9,10] have been identified using a 
combination of trichome-targeted EST analysis and 
biochemical assays. 

The precursors for bitter acid biosynthesis are 
branched-chain acyl-CoAs thioesters, which function as 
primers for polyketide synthesis, malonyl-CoA, which is 
the polyketide extender and dimethylallyl diphosphate 
(DMAPP), which serves as a prenyl donor. The acyl-CoAs, 
isovaleryl-CoA, isobutyryl-CoA and 2-methylbutyryl-CoA, 
are obtained by the degradation of the branched-chain 
amino acids (BCAA) leucine, valine and isoleucine, 
respectively (Figure la). The study of BCAA metabolism 
is of particular interest in hop since it is one of the few 
plants that produces large amounts of BCAA-derived 
natural products in specialized secretory structures. Some 
Solanaceae species utilize BCAAs in the trichome- 
localized production of acyl sugars [11,12]. 

BCAA biosynthesis in plants consists of eight enzymatic 
steps to produce leucine and valine from pyruvate, or 
isoleucine from threonine and pyruvate (Figure lb) 
(reviewed in Binder et al, [13]). Branched-chain amino- 
transferase (BCAT) enzymes catalyze both the final step 
in BCAA biosynthesis and the first step in BCAA 
degradation. These two phases are physically separated, 
with biosynthesis occurring in plastids and catabolism 
in mitochondria. BCATs have been studied in both 
Arabidopsis and tomato with each containing six and 
seven isoforms, respectively [14-16]. BCAA derived 
CoA-esters are generated by the branched-chain keto- 
acid dehydrogenase complex (BCKDH) which is made 
up of four subunits [17]. Three of these subunits, keto- 
acid dehydrogenase Ela and El|3, and dihydrolipoyl 
acyltransferase (E2), are specific to the BCKDH complex 
while the dihydrolipoyl dehydrogenase (E3) subunit also 
functions in other metabolic pathways [18]. 

Labeling studies have shown that the DMAPP 
required for the prenylation of bitter acids is provided 
by the methyl-D-erythritol 4-phosphate (MEP) pathway 
localized in plastids [19]. MEP-derived isoprenoids are 
also used in the biosynthesis of mono- and sesquiter- 
penes, which co-occur with bitter acids in lupulin 
glands [8]. The MEP pathway is initiated with 
glyceraldehyde-3-phosphate and pyruvate; the latter is 
also used in BCAA biosynthesis (Figure lb). Previous 
investigations have shown that transcripts for the seven 



MEP pathway enzymes and DMAPP/IPP isomerase are 
abundant in lupulin glands but their expression was not 
compared with leaves or other organs [5,8]. 

The final steps in bitter acid biosynthesis are unique 
to hop. The type III polyketide synthase enzyme valer- 
ophenone synthase (VPS) condenses a BCAA-derived 
acyl-CoA starter molecule with three molecules of 
malonyl-CoA to form the polyketide core (e.g. phlori- 
sovalerophenone, PIVP) [20] (Figure lb). One or more 
aromatic prenyltransferase enzymes perform two C5 
prenylations to yield diprenyl intermediates (e.g. 
deoxyhumulone) [21]. At this point the a- and |3-acid 
pathways diverge with a third prenylation forming the 
|3-acids while an oxygenation reaction, catalyzed by an 
as yet unknown humulone synthase enzyme, yields the 
a-acids [22]. 

Here we used RNA-seq to assemble and compare the 
transcriptomes of lupulin glands, cones (with glands 
removed) and leaves. We analyzed the expression of 
genes involved in precursor biosynthesis, including 
BCAAs, MEP-derived isoprenoids and malonyl-CoA, 
and the final steps in bitter acid formation. We focused 
on the BCAT enzymes which play key roles in the 
formation of branched-chain acyl-CoAs. Two BCAT- 
encoding genes were present, HlBCATl and HIBCAT2, 
with transcripts for both more abundant in lupulin glands 
compared with leaves and cones. HlBCATl and H1BCAT2 
were experimentally-localized using GFP fusion proteins 
to the mitochondrion and plastid, respectively. Assays 
of recombinant HlBCATl and H1BCAT2 showed they 
catalyze biosynthetic and catabolic reactions. These 
findings provide evidence for lupulin gland-specific 
BCAA and isoprenoid metabolism to produce precursors 
for bitter acid biosynthesis. 

Results 

Transcriptome sequencing and assembly 

Lupulin glands were isolated from the mature hop cones 
of four commercial cultivars grown for their high a-acid 
content: 'Taurus', 'Nugget', 'Magnum' and Apollo'. Leaves 
and cones were also collected from 'Taurus' and 'Apollo' 
cultivars. The cone samples consisted of the photosyn- 
thetic bracts and bracteoles that remained after the 
removal of lupulin glands. RNA purified from eleven 
samples (five gland, three cone, three leaf) was reverse- 
transcribed and the cDNAs sequenced on a single lane 
of an Illumina Hiseq 2000 instrument as a multiplexed 
sample. Each cDNA produced 15.7-24.7 million reads 
(Additional file 1: Table SI). The reads obtained from 
all samples were assembled using Trinity [23] and 
refined using Cdhit-EST [24] to produce a transcrip- 
tome of ~ 170 000 contigs with an average read length of 
745 bp. This transcriptome was used as the reference 
for RNA-seq analysis of the individual tissues using 



Clark et at. BMC Plant Biology 201 3, 1 3:1 2 
http://www.bionnedcentral.conn/1 471-2229/1 3/1 2 



Page 4 of 14 



CLC bio. At least 87% of the reads from each sample 
were mapped to the reference transcriptome for each 
tissue (Additional file 1: Table SI). 

Comparison of metabolic transcripts in lupulin glands, 
cones and leaves 

We compared the expression of mRNAs encoding 
primary metabolic enzymes (Calvin cycle, tricarboxylic 
acid (TCA) cycle and glycolysis) between the three 
tissues (Figure 2a, Additional file 2: Table S2). As 
expected, Calvin cycle genes were predominantly repre- 
sented in the photosynthetic tissue (leaves and cones) 
but not in the lupulin glands. Notable exceptions were 
the Calvin cycle enzymes involved in glycolysis, which 
are important for the flow of carbon from sucrose to 
pyruvate in the glands. An isoform of the Rubisco small 
subunit was one of the most abundant genes in the leaf 
transcriptome, while the Rubisco large subunit was not 
highly expressed. TCA cycle genes showed only minor 
differences in expression levels between the three 
tissues. Transcripts encoding glycolytic enzymes were 
also represented at comparable levels in leaves and 
lupulin glands, however, some genes within this 
pathway showed specificity to one tissue over the other. 
Several genes encoding enzymes required for sucrose 
utilization (e.g. fructokinase, neutral alkaline/ invertase 
and sucrose synthase) were more highly expressed in the 
gland versus leaf tissue. 

MEP pathway genes were more highly expressed in 
lupulin glands compared with green tissues. As 
expected from the need for malonyl-CoA for polyketide 
synthesis, transcripts encoding acetyl-CoA carboxylase 
[25] and ATP-citrate lyase [26] were abundant in the 
glands (Figure 2b). 

Most genes in the BCAA biosynthetic pathway 
showed higher expression in lupulin glands compared 
to cones and leaves (Figure 2b), with the transcript for 
ketol-acid reductoisomerase showing very high abun- 
dance. Exceptions were the low expression of threonine 
deaminase, which is only required for isoleucine biosyn- 
thesis, and the regulatory subunit of acetohydroxyacid 
synthase, which had minimal expression in all three 
tissues (Additional file 2: Table S2). Four putative 
BCAT-encoding genes were identified in the transcrip- 
tome. Two of these genes, denoted Humulus lupulus 
BCATl (HlBCATl) and HIBCAT2, were abundant in the 
lupulin glands and HlBCATl showed gland-specific 
expression (Figure 2b). Formation of the BCAA derived 
CoA-esters is carried out by the BCKDH complex 
(Figure lb). Two isoforms of the Ela subunit were iden- 
tified in the dataset and the expression of both isoforms 
along with the Elp and E2 subunits were increased in 
the lupulin glands (Figure 2b). The expression of the E3 
subunit was similar in all tissues. 



VPS was one of the most abundant mRNAs in the 
lupulin gland dataset, dwarfing all other genes involved 
in secondary metabolism examined here (Figure 2b). 
Two aromatic prenyltransferase genes, the previously 
characterized HlPTl [21] and HIPT2, showed similar 
gland-specific expression. 

Lupulin gland-specific transcription factors 

We also analyzed the RNA-seq data for transcription 
factors that may be involved in regulating bitter acid 
biosynthesis or other lupulin gland-specific metabolic 
activities. Forty-eight transcription factors were signi- 
ficantly more abundant in lupulin glands compared to 
the leaves (Additional file 3: Table S3). Of these, a subset 
of 29 was also higher in glands compared to cones. 
Eleven transcription factors were expressed at more 
than ten-fold higher levels in the lupulin glands com- 
pared to cones and leaves (Table 1). 

Validation of gene expression by quantitative RT-PCR 

The RNA-seq analysis of the BCKDH complex, lupulin 
gland-specific BCATs, VPS and HlPTl was confirmed 
using qRT-PCR. HlBCATl was lupulin gland specific 
while HIBCAT2 was expressed in both lupulin glands 
and leaves (Additional file 4: Figure SI). Ela, El^ and 
E2 of the BCKDH complex were more highly expressed 
in lupulin glands while the E3 subunit was represented 
equally in all tissues. The genes encoding VPS and 
HlPTl were highly expressed in the lupulin glands while 
essentially non-existent in the leaf tissue. All of these 
findings were in agreement with the RNA-seq data 
(Figure 2b). 

Metabolite profiling of BCAAs, acyl-CoAs and bitter acids 

We quantified BCAA levels in lupulin glands, cones and 
leaves using HPLC with fluorescence detection of their 
derivatized adducts. Leucine, valine and isoleucine 
levels were similar in lupulin glands and leaves, with 
valine the most prominent (Figure 3a). BCAA levels in 
the cone tissue were highly variable and, in the case of 
leucine and isoleucine, were also elevated compared to 
glands and leaves. 

Using a sensitive LC-MS method with multiple 
reaction monitoring, we measured the abundance of the 
acyl-CoAs that are used as substrates in bitter acid 
biosynthesis (malonyl-CoA, isovaleryl-CoA, isobutyryl- 
CoA and 2-methylbutyryl-CoA). Although acetyl-CoA 
is not used by VPS directly, it is converted to malonyl- 
CoA by acetyl-CoA carboxylase and therefore we included 
it in the analysis. Both acetyl-CoA and malonyl-CoA 
exhibit elevated levels in lupulin glands (Figure 3b). The 
starter substrates for the VPS reaction, isovaleryl-CoA, 
isobutyryl-CoA and 2-methylbutyryl-CoA, accumulated in 
glands while remaining largely absent in other tissues. 
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Figure 2 Scatter plot comparison of gene expression of metabolic pathways between lupulin glands, leaves and cones. Data is 
presented as reads per kb exon model per million mapped reads (RPKM). Corresponding gene names and expression values can be found in 
Table S2. a) Expression of genes encoding enzymes of primary metabolic pathways (Calvin cycle, TCA cycle and glycolysis), b) Expression of 
genes encoding enzymes of BCAA metabolism (BCAA, BCKDH complex, HIBCATl and HIBCAT2), MEP pathway, malonyl-CoA synthesis and bitter 
acid biosynthesis. 
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Table 1 Transcription factors that are ten-fold more 
abundant in the lupulin glands compared to leaves and 
cones 



Feature ID 




Description 


Fold change 








Gland^ vs 
Leaf*" 


Gland^ vs 
Cone'' 


comp76122_c0_ 


_seq4 


MYB 


498.2 


21.63 


comp50231_c0_ 


_seql 


WRKY 


138.5 


30.8 


comp68974_c0_ 


_seql 


WRKY 


127.2 


34.3 


comp53904_c0_ 


_seql 


MYB (CAN883 
homolog^) 


75.6 


14.5 


comp69082_c0_ 


_seql 


MYB 


67.3 


41.5 


connp70463_c0_ 


_seql 


AP2-like ethylene 
response 


56.6 


36.4 


comp58890_c0_ 


_seql 


Homeobox-leucine 
zipper 


54.9 


23.0 


connp27968_c0_ 


_seql 


Zinc finger 


50.0 


10.1 


connp72691_c0_ 


_seql 


Homeobox-leucine 
zipper 


45.3 


16.7 


connp73886_c0_ 


_seql 


WRKY 


39.9 


23.4 


connp63432_c0_ 


_seql 


MYB (CAN 738 
homolog^) 


18.8 


23.1 



Fold change represents the difference between mean RPKM values for gland 

compared to leaf and cone. 

^n=5. 

^n=3. 

' Marks et al. [38]. 



Isovaleryl-CoA was five times more abundant than the 
other branched-chain acyl-CoAs. 

HPLC analysis of the major bitter acids and their 
immediate precursors showed, as expected, that the 
a-acids accumulated to higher levels than the |3-acids 
and that humulone was the most abundant end-product 
(Figure 3c). These results correlated well with the levels 
of acyl-CoA precursors (Figure 3b). We also quantified 
the diprenylated precursors (e.g. deoxyhumulone) of the 
bitter acids. Small amounts of deoxyhumulone and 
deoxycohumulone were present while deoxyadhumu- 
lone could not be detected (Figure 3c, inset). The 
absence of deoxyadhumulone may reflect the overall 
lower levels of isoleucine-derived bitter acids. Although 
we attempted to identify the products from the VPS 
reaction such as phlorisovalerophenone, as well as 
the monoprenylated forms, these metabolites were 
undetectable. 

HIBCAT1 and HIBCAT2 catalyze aminotransferase 
reactions with BCAAs 

To further investigate BCAA metabolism in lupulin 
glands, we characterized the biochemical properties of 
the two gland-specific BCAT enzymes. HlBCATl is 393 
amino acids long with calculated molecular weight of 
43.2 kDa; H1BCAT2 is a 44.5 kDa (408 amino acid) 
protein. Nucleotide sequences for both enzymes are 



available in GenBank [HlBCATl, JQ063073; HIBCAT2, 
JQ063074]. The mature forms (i.e. lacking predicted 
transit peptides) of HlBCATl and H1BCAT2 were 
expressed as his-tagged recombinant proteins in E. coli 
and purified. Enzyme assays were performed in both the 
forward (anabolic, e.g. leucine to 2-ketoisocaproate) and 
reverse (catabolic, e.g. 2-ketoisocaproate to leucine) 
directions to test catalytic activities and to determine 
kinetic parameters (Table 2). Both enzymes functioned 
as BCATs and showed reversibility, as would be 
expected based on their ping pong reaction mechanism 
[27]. Unlike previous work in tomato, where the BCATs 
exhibited kinetic preferences corresponding to anabolic 
or catabolic functions [15], the hop BCATs had similar 
/<ni values in either direction. 

Subcellular localization of hop BCATs 

HlBCATl, H1BCAT2 and H1BCAT3 contain putative 
N-terminal signal peptides while H1BCAT4 does not. 
Target? [28] analysis predicted that HlBCATl and 
H1BCAT3 are localized to the mitochondria, and 
H1BCAT2 to the plastids. We tested the subcellular 
localization of the two BCATs with gland-specific 
expression, HlBCATl and H1BCAT2, using transient 
expression in Nicotiana benthamiana leaves. Hop leaves 
were not used for this experiment because they proved 
difficult to infiltrate with Agrobacterium solution. 
HlBCATl-GFP constructs with signal peptide (Figure 4) 
or signal peptide with mature protein (Additional file 5: 
Figure S2) co-localized with the mitochondrial MTRK 
marker [29]. The full-length H1BCAT2 failed to localize 
to a recognizable compartment, however its signal 
peptide fused to GFP colocalized with the plastid 
PTRK marker (Figure 4, Additional file 5: Figure S2). 
HlBCATl and H1BCAT2 constructs without signal 
peptide accumulated in the cytosol (Additional file 5: 
Figure S2). 

Phylogenetic analysis of hop BCATs 

A phylogenetic tree was constructed using BCAT protein 
sequences from hop, Arabidopsis, tomato. Cannabis 
sativa and other plants obtained using the Phytozome 
database [22]. We observed a separation of dicot BCATs 
into two clades, which could be classified as biosynthetic/ 
plastidial and catabolic/mitochondrial based on the 
presence of experimentally-localized proteins in each 
group (Figure 5). HlBCATl and 3 grouped into the 
mitochondrial BCAT clade and H1BCAT2 grouped into 
the plastidial clade, which agrees with their predicted 
localizations by TargetP and our analysis of intracellular 
localization (Figure 4). H1BCAT4 was grouped with the 
other BCATs which do not appear to localize to either 
the plastid or mitochondrion. Although the monocot 
BCATs clustered into two distinct groups that were 
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Figure 3 Metabolite analysis of lupulin glands, leaves and 
cones, a) Amounts of BCAAs (leucine, valine and isoleucine). 

b) Amounts of acyl-CoA precursors for bitter acid biosynthesis. 

c) Amounts of major bitter acid products. Inset shows amounts of 
diprenylated acylphloroglucinol intermediates. For all analyses, 
values represent mean ± standard error (n=5 for glands, n=3 for 
cone and leaf). 



Table 2 Kinetic parameters for recombinant HIBCAT1 and 
HIBCAT2 enzymes in anabolic (forward) and catabolic 
(reverse) directions 



HIBCAT1 
Anabolic: 

Glutamate 

Ketoisocaproate 

Ketoisomethylvalerate 

Ketoisovalerate 
Catabolic: 

2-Oxoglutarate 

Leucine 

Isoleucine 

Valine 
HIBCAT2 
Anabolic: 

Glutamate 

Ketoisocaproate 

Ketoisomethylvalerate 

Ketoisovalerate 
Catabolic: 

2-Oxoglutarate 

Leucine 

Isoleucine 

Valine 



Km (MM) ^ 

720 ± 55.59 
35 ± 3.27 
99 ± 10.75 
99 ± 6.53 

30 ± 5.39 
40 ± 9.20 
230 ± 26.94 
340 ± 90.94 



690 ±88.51 
120 ± 11.43 
160 ± 11.63 
200 ± 1 1 .72 

210 ±58.86 
280 ± 69.56 
250 ± 41.07 
130 ± 27.74 



Vmax [\xrc\o\lmqlm\r\) ' 

5.59 ± 0.13 
4.50 ± 0.04 
2.95 ± 0.09 
4.23 ± 0.08 

4.50 ± 0.16 
4.84 ± 0.23 
19.40 ± 0.85 
12.60 ± 0.96 



135.90 ±6.03 
119.20 ± 2.97 
84.68 ± 1 .70 
96.38 ± 1 .45 

36.92 ± 2.68 
43.43 ± 2.55 
88.08 ± 3.83 
23.29 ±1.16 



^mean ± standard error (n=3). 

positioned between the dicot clades, there is no experi- 
mental locaUzation data for monocot BCATs at this 
time. 

Discussion 

lllumina sequencing provides new genomic resources for 
hop 

Although the transcriptome of lupuUn glands has been 
extensively studied using ESTs obtained through Sanger 
sequencing [5,8,30] and several 454 pyrosequencing 
datasets are available in NCBI, an in-depth comparison 
of the lupulin gland transcriptome with green tissues 
(i.e. leaves, cones) has not been reported. Schilmiller 
et al [31] previously used 454 pyrosequencing to 
compare gene expression in tomato trichomes, leaves 
and stems, and RNA-seq analysis was used to identify 
terpene synthases in tomato trichomes [32]. Microarrays 
have also been used to compare the trichome and stem 
transcriptome of alfalfa [33]. The lllumina sequences 
reported here significantly expand the genomic resources 
available for hop, in particular leaves and cones for which 
limited sequences are available. The high number of 
contigs in the de novo transcriptome assembly probably 
relates to the presence of "broken contigs", as evidenced 
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Figure 4 Subcellular co-localization of BOAT signal peptides fused to GFP with organelle specific markers using transient expression in 
Nicotiana benthamiana leaves. The HIBCATl signal peptide fused to GFP colocalized witli tine mitocliondria! MTRK nnarl<er (upper panel) while 
the HIBCAT2 signal peptide-GFP fusion colocalized with the plastidial PTRK marker (lower panel). Plasmid constructs in Agro bacterium were 
agro-infiltrated into N. benthomiono leaves and fluorescent protein expression visualized using confocal microscopy. Scale bars represent 20 |jm. 



by the relatively short contig length (745 bp). This may 
have resulted from our use of combined RNA-seq data 
from several hop cultivars, the lack of sufficient sequen- 
cing depth and the presence of introns from alternative 
splicing events in some transcripts. In addition to the 
further development of transcriptome resources for 
hop, an important next step is to link the transcriptome 
with a whole genome sequence and genetic map. How- 
ever, the estimated 2.8 Gbp size of the Humulus lupulus 
genome [34] and the highly heterozygous nature of 
this dioecious species present obvious challenges for 
genome assembly. 

Lupulin glands synthesize precursor metabolites for bitter 
acid production 

The expression pattern for genes related to photosyn- 
thesis (e.g. Calvin cycle), and the absence of bitter acids 
from the cone minus gland samples (Figure 3c), indi- 
cates that separation of the non-photosynthetic glands 
from photosynthetic cones was largely complete. Ana- 
lysis of gene expression for primary metabolic pathways 
(TCA cycle and glycolysis) showed the reliance of lupu- 
lin glands on surrounding green tissues for the supply of 
sugars. The genes involved in the Calvin cycle were 
more abundant in the leaves and cones compared to 
lupulin glands (Figure 2a). The only exceptions to this 
pattern were genes encoding enzymes that function in 
both the Calvin cycle and glycolysis. The TCA cycle had 
similar expression levels across the leaf, cone and 



lupulin gland samples which can be attributed to its role 
in supplying intermediates for many pathways. Among 
the glycolysis transcripts three genes ifructokinase, 
neutral/ alkaline invertase and sucrose synthase) stood 
out for their high expression in the lupulin glands 
compared to leaves. All three are involved in sucrose 
mobilization, underscoring the exogenous supply of 
sugars to the lupulin gland. This is in agreement with 
the previous suggestion that sucrose would be an ideal 
energy source and precursor for secondary metabolism 
in glandular trichomes [6]. 

The high amounts of bitter acids produced by elite hop 
cultivars represents a large carbon commitment within 
the lupulin glands. Given that sugars are imported, we 
analyzed the RNA-seq dataset to determine if other 
precursors (e.g. BCAAs and prenyl diphosphates) are 
imported or synthesized in situ. Both the MEP and 
BCAA pathways were more abundantly represented in 
the lupulin glands (Figure 2b), a finding that provides 
evidence for the de novo biosynthesis of BCAAs and 
DMAPP in these cells. BCAA biosynthesis in lupulin 
glands is logical given the demand for the acyl side- 
chains of the BCAAs but not the amino group, and the 
need to prevent the accumulation of excess nitrogen. 
This finding is in agreement with Schilmiller et al. [6] 
who noted that while phenylpropanoid-derived metabo- 
lites are synthesized in trichomes (e.g. phenylpropenes 
in basil) the import of phenylalanine into glandular 
trichomes was minimal. 



Clark et at. BMC Plant Biology 201 3, 1 3:1 2 
http://www.bionnedcentral.conn/1 471-2229/1 3/1 2 



Page 9 of 14 




Figure 5 A phylogenetic tree of BCAT proteins in plants constructed using the neighbour-joining method. BCATs group into plastidial 
(orange line) and mitocliondrial (blue line) clades, with apparent division into dicot and monocot subclades. A third clade (grey line) includes 
BCATs that fall into neither group. BCATs with experimentally-confirmed organelle localization are labelled: P, plastid; M, mitochondrion; C, 
cytosol. Species abbreviations are found in Methods (At, Arobidopsis tholiono; HI, Humulus lupulus; SI, Solonum lycopersicon). *As we note in the 
Discussion, the plastid localization of AtBCAT2 requires further analysis. 



Transcripts corresponding to the enzymes of the 
polyketide portion of the bitter acid pathway (e.g. VPS) 
were among the most abundant in the glandular trich- 
ome dataset (Figure 2b). The presence of two highly 
expressed aromatic prenyltransferases is notable since 
HlPTl has only been shown to catalyze the first preny- 
lation step in bitter acid biosynthesis [21]. The a- and 
|3-acids require two and three prenyl groups, respect- 
ively, and the second prenyltransferase enzyme 
(H1PT2, 45% amino acid identity to HlPTl) may 
catalyze these steps. 



Lupulin gland-specific expression of transcription factors 

RNA-seq analysis also revealed a number of lupulin 
gland-specific transcription factors which may be 
involved in regulating bitter acid biosynthesis or other 
metabolic or developmental processes in glandular 
trichomes (Table 2). Matoucek and co-workers have 
identified several transcription factors that may control 
flavonoid biosynthesis in hop through the regulation of 
chalcone synthase (CHS) (e.g. HlMYBl, H1MYB2, 
H1MYB3, H1MYB7, HlbHLH2 and HlWDRl [35-37]). 
Although all of these transcription factors could be 
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identified in the hop transcriptome reported here, only 
transcripts of HIMYB7 and HlbHLH2 were more abun- 
dant in lupuUn glands compared to leaves (Additional 
file 3: Table S3). This expression pattern makes sense 
since CHS, which plays a role in both the gland-specific 
production of xanthohumol and in the formation of 
flavonoids, is expressed throughout the plant. Indeed 
our RNA-seq analysis shows that CHS_H1 is expressed 
in lupulin glands and leaves (Additional file 2: Table S2), 
with the three-fold higher expression levels in glands 
likely reflecting the production of xanthohumol. Two 
MYB transcription factors, CAN883 and CAN778, have 
been shown to be highly expressed in cannabis glandular 
trichomes and implicated in regulating trichome-specific 
metabolic processes [38]. Both have close homologs to 
hop transcription factors that are over ten-fold more 
abundant in the lupulin glands (Table 2). CAN738, 
which has 70% amino acid identity to the hop MYB 
protein encoded by comp63432_c0_seql, was proposed 
to be involved in regulating cannabinoid formation. 
Since Cannabis and Humulus are closely related genera 
with glandular trichome-localized polyketide pathways, 
comp63432_c0_seql may be a candidate for regulating 
the bitter acid pathway. 

Lupulin glands contain both mitochondrial and plastidial 
BCATs 

The biosynthesis and subsequent degradation of the 
BCAAs requires both metabolic and catabolic BCATs. 
The lupulin gland datasets contained two BCATs, 
HlBCATl and HIBCAT2, which were abundant in the 
glandular trichomes. Based on RNA-seq and qRT-PCR 
analysis (Figure 2b, Additional file 4: Figure SI) 
HlBCATl is exclusively expressed in the lupulin glands. 
Transient expression of the HlBCATl-GFP fusion con- 
struct confirmed its mitochondrial targeting (Figure 4). 
H1BCAT2 with N-terminal signal peptide failed to 
localize to the plastid and appeared to form aggregates 
in the N. benthamiana epidermal cells (Additional file 
5: Figure S2), however, the predicted H1BCAT2 signal 
peptide fused to GFP was targeted to the plastid 
(Figure 4). These locaUzations were further supported 
through the construction of a phylogenetic tree compar- 
ing BCATs from a variety of plant species (Figure 5). 
Problems with the experimental localization of hop 
plastid proteins have been previously reported for a 
GPP synthase small subunit (GPP-ssu) identified in 
lupulin glands [8]. 

Kinetic analysis suggests lack of specialization for hop 
BCATs 

Kinetic analysis has provided equivocal evidence for the 
catalytic specialization of BCATs for BCAA biosynthesis 
or catabolism. In tomato, biosynthetic BCATs (S1BCAT3 



and 4) had values with the keto acids that were an 
order-of-magnitude less than the catabolic enzymes 
(SlBCATl and 2), indicating specialization of the former 
for forward (biosynthetic) reactions [15]. In contrast, 
Arabidopsis AtBCATl, which has been shown to be a 
catabolic enzyme, had values for the keto acids 
(0.036-0.8 mM) that were much lower than for BCAA 
substrates (0.4-6 mM), suggesting a preference for the 
biosynthetic reaction. Indeed, the values for AtBCATl 
with the keto acids were lower than those reported with 
the plastidial, biosynthetic enzyme AtBCAT3 [39]. The 
kinetic properties of the hop BCATs appear to be more 
similar to the Arabidopsis enzymes. The catabolic 
enzyme HlBCATl exhibited lower values for keto 
acids, and therefore an apparent preference for the for- 
ward (biosynthetic) reaction, compared with than 
H1BCAT2 (Table 2). It is noteworthy that HlBCATl 
showed a lower for leucine compared with valine 
and isoleucine which may reflect this enzyme s preferred 
role in synthesizing isovaleryl-CoA for humulone 
biosynthesis. The lower V^ax values for the HlBCATl 
may also explain why it is more highly expressed 
compared to H1BCAT2 in the lupulin gland tissue 
(Figure 2). 

Phylogenetic analysis indicates separation of biosynthetic 
and catabolic BCATs 

To further analyze the function of the hop BCATs, we 
constructed a phylogenetic tree of diverse BCATs in- 
cluding well characterized enzymes from Arabidopsis 
and tomato. This analysis was successful in separating 
biosynthetic (plastidial) and catabolic (mitochondrial) 
BCATs into separate clades. BCATs are also reported to 
accumulate in alternative cellular compartments such as 
the cytosol and vacuole for an as of yet unknown func- 
tion [15,16] and these genes likely make up the third 
group that includes several cytosolic enzymes and 
H1BCAT4 (Figure 5). An anomalous grouping was 
observed for AtBCAT2, which has been shown to be 
plastid localized [6] but which was placed in the mito- 
chondrial clade. A previous study found AtBCAT2 to be 
the only Arabidopsis gene to be co-expressed with other 
genes involved in BCAA catabolism [40], a role one 
would expect to be occupied by a mitochondrial BCAT. 
In addition, TargetP analysis of AtBCAT2 predicts mito- 
chondrial rather than plastidial localization. 

Conclusions 

This study provides an in-depth analysis of the tran- 
scriptional activity of primary and specialized metabol- 
ism leading to hop bitter acids. RNA-seq analysis 
showed that both the MEP and BCAA pathways were 
highly expressed in lupulin glands together with genes 
required for remobilization of carbon from sucrose. 
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Branched-chain acyl-CoAs and bitter acids are present 
in higher levels in isolated lupulin glands compared 
with cones (with glands removed) and leaves. A plasti- 
dial BCAT enzyme involved in BCAA biosynthesis and 
a mitochondrial BCAT catalyzing BCAA degradation 
both show lupulin gland-specific expression. Together 
these findings indicate that the glandular trichomes are 
the site of biosynthesis for the precursors for bitter 
acid biosynthesis. The deep transcriptome sequencing 
reported here significantly expands the genomic 
resources available for hop. This new information, and 
our analysis of BCAA metaboUsm, will be useful for 
the further elucidation of enzymes and regulatory pro- 
teins involved in bitter acid production and for breed- 
ing hop cultivars with increased or modified bitter acid 
content. 

Methods 

Hop transcriptome assembly and RNA-seq analysis 

RNA was isolated from lupulin glands sampled from 
hop cultivars 'Taurus', 'Nugget' and 'Magnum' (collected 
in Mainburg, Germany in August 2005), from 'Taurus' 
(collected in Saskatoon, SK in September 2011) and 
'Apollo' (collected in Yakima, WA in August 2009). 
RNA was isolated from leaves and cones-minus-glands 
from 'Apollo' (collected Yakima 2009) and 'Taurus' 
(collected Mainburg 2005 and Saskatoon 2011). Glands 
were separated from cones by stirring in liquid nitrogen 
followed by filtration through a 1 mm screen to separate 
cone material. Total RNA was isolated using a CTAB 
method and an RNeasy kit with on-column DNasel 
treatment (Qiagen) [41]. RNA was prepared for sequen- 
cing by the National Research Council Canada DNA 
Technologies Unit using the lUumina TrueSeq RNA 
sample preparation platform v.2 with multiplex labeling 
and subjected to paired-end sequencing using an 
lUumina HiSeq 2000 with a read length of 101 nt. 
Trimmomatic 0.15 [42] was used to remove adaptor 
sequences and trim bases with quality lower than 20 
(Phred 33 quality scores). The transcriptome assembly 
was performed with Trinity [23] and refined using cd- 
hit-EST [24] to combine identical reads. CLC bio v.5 
was used for RNA-seq analysis and statistical analysis 
via the Baggerly's test [43]. Annotation of the transcrip- 
tome was performed using Blast2GO [44]. 

The Illumina data for all samples have been deposited in 
NCBI under Bioproject PRJNA175602 with the transcrip- 
tome available via accession number GAAWOOOOOOOO. 
RPKM values for all transcripts is found in Additional file 
6: Table S4. 

qRT-PCR analysis 

cDNA was synthesized from 0.5 (ig total RNA from 
lupuUn glands of 'Taurus' and 'Nugget' sampled in 



Mainburg, Germany and 'Taurus' sampled in Saskatoon, 
Canada. RNA was isolated from leaves from 'Taurus' 
sampled in Mainburg, Germany in 2005, Saskatoon, 
Canada and 'Apollo' from Yakima, USA using the 
QuantiTect reverse transcription kit (Qiagen). Primers 
(Additional file 7: Table S5) were designed to amplify 
products of 101-110 bp. PCR reactions were performed 
with a StepOnePlus real-time PCR system (Applied 
Biosystems) using the QuantiFast SYBR Green qPCR 
kit (Qiagen). |3-tubulin was used a reference gene. 
StepOnePlus software was used to calculate AACt 
values [45]. 

Metabolite analysis 

For BCAA analysis, samples were weighed and ground 
in 5 vol (based on 1 mL per mg fresh weight of plant 
material) of 0.1 N HCl, centrifuged at 20,000 g for 10 
min and derivatized using the AccQ-Tag reagent 
(Waters). Derivatized amino acids were analyzed with 
a Waters 2695 separation module using a Waters 
Symmetry C18 column (5 (im, 2.1 x 100 mm) and a 
Waters 2475 fluorescence detector as described by 
Cohen and Michaud [46]. Amino acids were quantified 
using standard curves. 

For CoA analysis, samples were ground in liquid 
nitrogen and extracted twice with 10 vol of 5% TFA with 
internal standards (^^Cs malonyl-CoA, butenoyl-CoA 
and benzoyl-CoA). Samples were evaporated, resus- 
pended in 900 (iL of acetonitrile: isopropyl alcohol: 100 
mM KH2PO4: aqueous acetic acid (9:3:4:4, v/v/v/v) and 
solubilized by vortexing and sonication. Oasis WAX 
weak anion exchange cartridges (Waters) were primed 
with methanol, equilibrated with acetonitrile: isopropyl 
alcohol: water: acetic acid (9:3:4:4, v/v/v/v), the samples 
loaded and washed with the equilibration mixture. 
Samples were eluted with 250 mM ammonium hydroxide 
in 80% methanol, dried and stored at -20°C. The resi- 
due was reconstituted in 100 (iL of water: acetonitrile 
(95:5, v/v) containing 5 mM triethylamine (TEA) and 
3 mM acetic acid, vortexed, sonicated and centrifuged. 
The injection volume was 10 (iL. CoAs were separated 
on a Waters Acquity UPLC system with a Waters BEH 
CIS column (1.8 (im, 2.1 x 100 mm) with a Vanguard 
guard column. Mass detection used a Quattro Ultima 
(Waters) in positive ion, multiple reaction monitoring 
(MRM) mode. Solvents were A=5 mM TEA and 3 mM 
acetic acid in water, and B=5 mM TEA and 0.3 mM 
acetic acid in acetonitrile: water (95:5, v/v). A flow rate 
of 200 (iL/min was used with an initial solvent compo- 
sition of 98:2 A:B changing to 87:13 A:B over 19 min. 
Detection of the acyl-CoAs using MRM involved 
neutral loss of 507 Da from each protonated molecular 
ion [47]. Quantification by standard curves used 
MassLynx 4.1 (Waters). Since there was no standard for 
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2-methybutyryl-CoA, quantification was based on the 
calibration curve constructed for n-valeryl-CoA. 

For polyketide analysis, samples were ground in 10 vol 
of 80% MeOH, filtered through a 0.45 (im Spin-X filter 
(Corning) and separated on a Waters 2695 separation 
module using a Waters Symmetry C18 column (5 (im, 
2.1 X 100 mm). A flow rate of 0.3 mL/min was used with 
an initial solvent composition of acetonitrile: water 
(46:54, v/v) with 0.05% formic acid changing to aceto- 
nitrile: water (95.5:4.5, v/v) with 0.05% formic acid over 
25 min. Metabolites were detected using a Waters 
photodiode array detector at 237 nm and 3100 Mass De- 
tector in SIR with positive electrospray and a cone volt- 
age of 30 V. Metabolite levels were determined with 
standard curves of deoxyhumulone, humulone and 
lupulone. 

BCAT characterization 

Oligonucleotides used in cloning procedures are found 
in Additional file 7: Table S5. Signal peptides were 
predicted using Mitoprot [48] and ChloroP [49]. The 
open-reading frames of HlBCATl and H1BCAT2 lacking 
signal peptides were amplified using Pfu Ultra II poly- 
merase (Agilent), the amplification products cloned into 
pCR8/GW/TOPO (Invitrogen) and recombined into 
pHIS8/GW using LR Clonase (Invitrogen). HlBCATl 
and H1BCAT2 in pHIS8/GW were transformed into E. 
coll BL21(DE3) Rosetta2 pLysS cells (Novagen). Cells 
were grown in Luria-Bertani media containing 50 [ig/L 
kanamycin and 34 (ig/L chloramphenicol to an ODgoo of 
0.6 at 37°C, cooled to 18°C, induced with 0.3 mM IPTG 
and incubated for 4 h at 18°C. Bacterial cultures were 
pelleted and frozen at -80°C. Pellets were thawed, resus- 
pended in 20 mM Tris HCl (pH 8.0), 300 mM NaCl, 100 
[iM pyridoxal-5-phosphate (PLP), 10% glycerol, 5 mM 
imidazole, 5 mM p-mercaptoethanol, 0.2% triton X-100, 
0.1% CHAPS, 0.2% sarkosyl, 0.5 mM PMSF, 1 (ig/ml 
pepstatin A and 1 (ig/ml leupeptin, and lysed by sonic- 
ation. Recombinant proteins were purified from the cell 
lysate using Talon resin (Clontech). The resin-bound 
protein was washed with 20 mM Tris HCl (pH 8.0), 300 
mM NaCl, 10% glycerol, 10 mM imidazole, 5 mM p- 
mercaptoethanol and eluted with the same buffer com- 
position containing 150 mM imidazole. Samples were 
desalted and concentrated using Amicon Ultra filters 
(Millipore) and stored at -80°C. Protein concentration 
was determined using the Bradford assay [50] and pro- 
tein purity was assessed at >95% via SDS-PAGE gel ana- 
lysis with Coomassie brilliant blue stain. 

Anabolic reactions were linked to glutamate dehydro- 
genase to allow monitoring with NADH turnover as 
described by Prohl et al [51]. Anabolic assays contained 
66 (ig/ml HlBCATl or 12 (ig/ml H1BCAT2, 250 mM Tris 
HCl (pH 8.5), 50 (iM PLP, 1.5 mM dithiothreitol, 0.2 mM 



EDTA, 100 mM ammonium chloride, 0.2 mM NADH, 3 
units/mL glutamate dehydrogenase (Sigma), 2 mM glu- 
tamate for branched-chain keto-acid kinetics and 0.5 mM 
ketoisovalerate for glutamate kinetics. Assays were run in 
triplicate at 22°C and observed in real time at 340 nm with 
an Ultraspec 3100 pro spectrophotometer (Biochrom). 
Time point assays were employed for catabolic reactions. 
Catabolic assays contained 5.4 (ig HlBCATl or 1.2 (ig 
H1BCAT2, 20 mM Tris HCl (pH 8.5), 50 (iM PLP, 1.5 mM 
dithiothreitol, 0.2 mM EDTA, 2 mM 2-oxoglutarate for 
BCAA kinetics and 2 mM leucine for 2-oxoglutarate 
kinetics. Assays were run in triplicate at 22°C for 2 min 
and halted with the addition of HCl. Samples were centri- 
fuged at 20,000 ^ for 10 min, derivatized using the AccQ- 
Tag reagent (Waters) and analyzed as described above. 

Analysis of subcellular localization 

Three constructs in pCR8/GW/TOPO were prepared for 
each gene: full-length including predicted signal peptide, 
the mature protein lacking the signal peptide, and the signal 
peptide only. Each was recombined into pEarleyGatel03 
to generate C-terminal GPP fusion constructs under the 
control of the 35S promoter [52]. Plasmid constructs of 
HlBCATs, and mitochondrial (MTRK) or plastid (PTRK) 
markers [29], were transformed into Agrobacterium 
C58 and infiltrated into Nicotiana benthamiana leaves as 
previously described [53]. Plants were grown for 2-4 days 
and the leaves directly observed with a Carl Zeiss Axiovert 
200M/ConfoCor 2/LSM 510 confocal laser scanning 
microscope. GFP was observed with an argon laser with 
excitation at 488 nm and emission at 505-530 nm. 
mCherry was observed with a helium laser with excitation 
at 543 nm and emission at 560-615 nm. 

Phylogenetic tree 

BCAT protein sequences for were obtained from the 
Phytozome database [22] for Manihot esculenta, 
Populus trichocarpa, Cucumis sativus, Prunus persica, 
Medicago truncatula, Glycine max, Carica papaya, 
Citrus Clementina, Eucalyptus grandis, Mimulus guttatus, 
Sorghum bicolor, Zea mays, Setaria italica, Oryza 
sativa, and Brachypodium distachyon. Previously char- 
acterized Arabidopsis and tomato proteins were 
included along with those from the Cannabis sativa 
genome [54]. All sequences were aligned using ClustalW 
[55] with variable N-terminal regions removed when 
present. Mega5 [56] was used to construct a neighbour- 
joining tree with 10,000 bootstrap replicates. 

Additional files 



Additional file 1: Table SI. Total lllumina sequencing reads for each 
RNA sample and the number of reads that mapped to the reference 
transcriptome. 
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Additional file 2: Table S2. RNA-seq analysis of metabolic genes 
compared between lupulin gland and leaf, lupulin gland and cone 
(minus glands), and leaf and cone (minus glands). Data is represented in 
reads per kb exon model per million mapped reads (RPKM). Note that 
each tissue comparison is located on a separate Excel sheet. 

Additional file 3: Table S3. Transcript abundance of transcription 
factors compared between lupulin gland and leaf, lupulin gland and 
cone (minus glands) and leaf and cone (minus glands). Data is 
represented in RPKM. 

Additional file 4: Figure SI. Relative abundance of select metabolic 
transcripts in lupulin gland and leaf samples as determined by qRT-PCR. 
Data is presented using relative quantitation (RQ) based on AACj values. 
Samples are normalized using (3-tubulin. Data represents the mean and 
standard error, n=3. Transcript abundance between panels a, b and c is 
not directly comparable in this figure, a) qRT-PCR analysis of HIBCATl and 
HIBCAT2. A single gland sample was used as reference, b) qRT-PCR 
analysis of genes encoding enzymes of the BCKDH complex (ketoacid 
dehydrogenase El a and Eip, dihydrolipoyi acyltransferase (E2), and 
dihydrolipoyi dehydrogenase (E3)). A single leaf sample was used as 
reference, c) qRT-PCR analysis of VPS and HIPTl. A single leaf sample was 
used as reference. 

Additional file 5: Figure S2. Subcellular co-localization of GFP fusion 
proteins of HIBCATl and HIBCAT2, with and without N-terminal signal 
peptides, compared with organelle specific marker proteins using 
transient expression in Nicotiona benthomiana leaves. MTRK and PTRK are 
mitochondrial and plastidial makers. Scale bars represent 20 [im. 

Additional file 6: Table S4. RPKM values for lupulin glands, cones 
(minus glands) and leaves. 

Additional file 7: Table S5. Sequences of oligonucleotides used. 
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